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Events vs Flux/Dose 

• NASA High-Z and E Transport (HZETRN) code calculates the 
average flux and dose of particles behind spacecraft and tissue 
shielding. 

• Monte-Carlo transport codes (GEANT, FLUKA, etc.) are 
cumbersome and not used for biophysics applications. 

• An event refers to the correlated energy depositions in time and 
space of cosmic ray interactions with cells controlled by the tissue 
matrix (environment). 

- Time-dependent transport codes are needed due to cell & 
tissue signaling activation and relaxation times: 

• Biological steady-state is altered by proton hits pre-, 
during, or post- HZE events. 

- Transport code must describe temporal and micro-spatial 
density of functions to correlate DNA and oxidative damage 
with non-targeted effects (signals, bystander, or other). 
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Biological Process Relaxation Time 

- Multiple events by GCR and SPEs for given process - 



DNA repair/ATM Relaxation Times 3-24 hrs 



Cucinotta et al. , Biochemical Kinetics of DSB Repair and Induction of y-H2AX Foci by Non-homologous End Joining, 
Rad. Res., 169, 214-222, 2008 


4 




Biological Process Relaxation Time 

- Multiple events by GCR and SPEs for given process - 

TGFp-SMAD Relaxation Times 1-5 days 

Figure 1. Schematic representation of Smad dependent TGF-p signaling pathway. 
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PLoS ONE 2(9): e936. doi:10.1371/journal.pone.0000936 
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Biological Process Relaxation Time 

- Multiple events by GCR and SPEs for given process - 


Repopulation Relaxation Times 

Differentiation 

Senescence 



1-30 days 


Cucinotta and Dicello, On the Development of Biophysical Models for Space Radiation Risk Assessment, NTRS 
20000083882 , 1999 


Transport Codes for Stochastic Models of 

Radiation Risks 



• New approaches to risk assessment will require event based 
models of particle transport that track time and spatial 
dependent interactions of particles in tissue structures 

• The GCR Event Based Risk Model (GERMcode) is a Monte-Carlo 
based approach for this purpose that builds on the success of 
HZETRN/BRYNTRN codes using QMSFRG 

• The GERMCode will incorporate stochastic distribution of 
incident particles 

v' Bi-directional transport allows to use FISHbowl spacecraft and organ 
geometry ray tracing 

■S Angular corrections can be added for small tissue samples where risk 
models are formulated using the stochastic approach 

■S The GERMcode will tally time-dependent events in support of new 
approaches to biological response models 




Heavy Ion Reactions 


Abrasion=projectile-target overlap 
(n, p, and cluster knock-out) 
Ablation = pre-fragment decay 
(n, p, d, t, h, alphas de-excitation) 
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NSRL for Biophysics Applications 


Approximate Composition 

^ 101 . 7 ^ 33 . 1^^36 

Density: 0.00194 g/cm 3 
Thickness: 1.2166 g/cm 2 
N: 2.09X10 22 atoms/g 
O: 6.81X10 21 atoms/g 
Al: 7.41X10 21 atoms/g 
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NSRL Bragg Curve Comparison to 
GCR Event-based Risk Model (GERMcode) 


56 Fe(0.59 GeV/u) 



48 Ti (0.98 GeV/u) 



Cucinotta FAetal., Radiation Protection Dosimetry, 143, 384-390, 2011, 
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NSRL Bragg Curve Comparison to 
GCR Event-based Risk Model (GERMcode) 
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Cucinotta FAetal., Radiation Protection Dosimetry, 143, 384-390, 2011 
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Thick Target Comparison with NASA's 

GERMcode 

Iron (1 GeV/u) on Polyethylene 



Cucinotta FAetal., Radiation Protection Dosimetry, 143, 384-390, 2011 
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Summary of GERMcode Accuracy for Physics 

• Atomic variables agree with experiments to within ±5% 

■S LET, Range, Straggling parameters 

• Absorption X-section within ± 5% 

• Elemental fragment X-section within ± 25% to H.l. experiment 

■S Errors are local in Z and E minimizing their impact 

• Comparison to NSRL Data: Excellent agreement at all depths 

• QMSFRG X-sections in HZETRN/GERMcode; classical X-sections 
in GEANT4, FLUKA and PHITS models 

• Focus of future work to add mesons and photons/electrons 
and to improve accuracy of event generator for light particles 

• Systems biology approach to risk prediction requires event 
based physical/biological models to account for stochastic 
transition rates 
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Future Plans 



• GCR Event-based Risk Model (GERM) wraps the existing physics 
code developed on the span of decades into a user-friendly 
graphics interface based on a fast Monte-Carlo algorithm 

• Radiation transport in GERMcode is based on Monte Carlo 
method to solve the transport problem for a distribution of 
particles present in space and track the evolution of individual 
particles within a material 

- The Monte Carlo approach will work together with the bi-directional ray-tracing 
technique following approach of HZETRN/BRYNTRN codes 

• From radiation transport in matter we advance the model to 
the cell and tissue effects that can address risk models that go 
beyond dose and dose equivalent 

• New technology (GPU) will enable the model to address full 
GCR simulations 
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NSRL GERMcode GUI vl.l 
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NSRL GERMcode GUI Overview 


• A stochastic simulation tool using track structure and nuclear 
interactions provides the description and integration of physical 
and biophysical events from mono-energetic ions. 

• A stochastic Monte-Carlo based model of radiation transport in 
spacecraft shielding and tissue is developed with the quantum 
multiple scattering model of heavy ion fragmentation (QMSFRG) 
and the energy loss processes. 

• For the scientists who participate in NSRL experiments or in data 
interpretation of such experiments, GERMcode provides the 
ability to: 

v' Model the beam line, shielding of samples and sample holders 

S Estimate basic physical and biological outputs of the designed 
experiments 
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GUI for the NSRL GERMcode 
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Radiobiological Model 


Radiobio property output: 
- Biological damage 
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Mono-Energetic Beam 

( Start ) 


A, Z, E, Material, 
Dose, 

Cell area, 

Radiobiological Model 
DNA volume 








Physical/radiobio 
property outputs: 

- LET Range curve 

- Nuclear extinction 

- Radial dose 

- Biological damage 


Nuclear interaction 
property output: 

- Probability of hits 


Track structure 
property output: 

- Energy deposit in 
DNA ' 


ft 




ft 





< End ) 
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Mono-Energetic Beam Transport 

( Start ) 
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Beam 
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Radiobiological Model 


Radiobio property output: 
- Biological damage 
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Parameter 

Charge number, Z 
Mass number, A 
Beam energy, 
MeV/u 

Material 

Dose, Gy 
Cell area, pm 2 
DNA volume 
[d x / of cylinder 
volume, unit in nm ) 


Transport Depth 


Radiobiological 

model 


User Input Control Parameters 


Mono-energetic beam 


Radiation transport in 
thick target 

1-28 

1-28 

1-58 

1-58 

50 - 1500 MeV/u 

50 - 1500 MeV/u 

Water 

Aluminum 

Polyethylene 

C0 2 

Graphite Carbon 

Water 

Aluminum 

Polyethylene 

C0 2 

Graphite Carbon 

0.0 -5.0 Gy 


0.1-1000 pm 2 



DNA segment (2 x 2) 
Nucleosome (10 x 5) 
DNA fiber (25x25) 




• Fixed Depth 

• Bragg Curve — 

• Biological Models 

• No Katz model 

• Cell survival 

• Chrom. aberration 

• Cell mutation 

• Mouse tumor model 

• No radiobio model 

• Cell survival 

• Chrom. aberration 

• Cell mutation 

• Mouse tumor model 






• Beam transport for fixed depth 


• Beam transport for Bragg curve depth 



Mouse: 



Longitudinal placement along the beam 



Mouse: 



Transverse placement to the beam 



Rat: 


• Beam 

Longitudinal placement along the beam 


Rat: 


transport 

Transverse placement to the beam 


in bio- 

Ferret : 


logical 

Longitudinal placement along the beam 


Ferret : 


sample 

Transverse placement along the beam 



T-25 flask 



T-75 flask 



Flaskette 



Chamber slide 



6-well plate 
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Ion Types in GERMcode with Default Nuclei Highlighted 
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Physical and Biophysical Properties of 
Mono-energetic Beams 


Input Parameter 


Mono-energetic beam 

Charge number , Z 

1-28 

Mass number, A 

1-58 

Beam energy, E(MeV/u) 

50 - 1500 MeV/u 

Material 

Water 

Aluminum 

Polyethylene 

C0 2 

Graphite Carbon 

Dose, Gy 


0.0 -5.0 Gy 

Cell area, p,m 2 


0.1-1000 )u,m 2 

DNA volume 

[d x / of cylinder volume, 

unit in nm) 


DNA segment (2 x 2) 
Nucleosome (10 x 5) 
DNA fiber (25 x 25) 


Radiobiological model 
(Katz model 1 ) 


• No Katz model 

• Cell survival 

• Chrom. aberration 

• Cell mutation 

• Mouse tumor model 


Output Input Parameter 


Physical property 

LET and Range curve 

(Z, A, E) Materia | 

Nuclear extinction 

(Z, A, E) Material 

Radiobiological 

property 

Radial dose 

(Z, A, E)jj ssue 

Biological damage of 
Katz model 1 

(Z, A, E) Tissue 

Nuclear 

interaction 

property 

Probability of hits 2 

(Dose, Cell Area, Z, A, E) Tissue 

Track structure 
property 

Energy deposition in 
DNA volume 3 

(Z, A, E) Tissue 


1 Biological damage using Katz model: 

Cell survival; 

Chromosomal aberration; 

Cell mutation; 

Mouse tumor model 

Probability of hits for Poisson distribution of ion hits 
3 DNA volume of d x / cylinder volume in nm\ 

DNA segment (2 x 2); 

Nucleosome (10x5 for 160 BP); 

Chromosome fiber (25 x 25) 
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Biophysical and Radiobiology Properties of NSRL Beam 

Transport 

(Monte-Carlo trials along path of primary ion) 


Input Parameter 

Charge number, Z 
Mass number, A 
Beam energy, E(MeV/u ) 

Material 


Transport depth 


Radiobiological model 


Radiation transport in 
thick target 


1-28 

1-58 

50 - 1500 MeV/u 
Water 
Aluminum 
Polyethylene 
C0 2 

Graphite Carbon 


• Fixed depth 

• Bragg curve depth 

• Biological models 


No Katz model 



Bio- 

physical 

property 


Radio- 

biological 

property 



Input parameter 

Output 

Fixed depth 

Bragg curve 
depth 

Biology 

model 

depth 

Depth-dose 

(Z,A,E) Material 

(Z,A,E) Material 

(Z, A, E ) T j Ssue 

Charge 

distribution 

(Z,A,E) Material 

(Z,A,E) Material 

(Z,A,E)j issue 

Multiplicity 
of events of 
HI, Neutron, 
Proton, or a 

(Z,A,E) Material 

(Z,A,E) Material 

(Z,A,E) Tissue 

Biological 
damage of 
Katz model * 1 

(Z,A,E)j jssue 

(Z,A,E)j jssue 

(Z,A,E)j jssue 


Cell survival 
Chrom. aberration 
Cell mutation 
Mouse tumor model 


1 Biological damage using Katz model: 
Cell survival; 

Chromosomal aberration; 

Cell mutation; 

Mouse tumor model 
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From Radiation Transport in Materials 
To Radiation Effects in Astronauts 


GERMcode results are applicable to biological events on 
the cell/tissue level. 

Energy imparted from a particle at certain material depth 
can be scored per pixel within a cell and per cell in a 
tissue matrix. 

Radiation transport in matter is applied to study tissue 
radiation effects within a human body. 

The scored stochastic biological effects can be DNA 
double strand breaks, apoptotic cells, or other processes. 

To speed up Monte Carlo simulations, a new technology 
refereed as General-Purpose Graphic Processor Unit 
(GPGPU) will be implemented. 


Homework 
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Beam Ion Species and Energies Used Previously at NSRL 


Ion Species 

Energy, 

MeV/u 

Maximum Intensity, 
ions per spill 

LET for a Water, 
KeV/pm 


50 - 2500 

6.4 x 10 11 

1.26-0.21 

4 He 

50 - 1000 

0.88 x 10 10 

5.01-0.89 

12 C 

65 - 1000 

1.2 x 10 10 

36.79-8.01 

16 Q 

50 - 1000 

0.4 x 10 10 

80.5-14.24 

20 Ne 

70 - 1000 

0.1 x 10 10 

96.42- 22.25 

28 Si 

93 - 1000 

0.3 x 10 10 

151-44 

35 CI 

500 - 1000 

0.2 x 10 10 

80-64 

40 Ar 

350 

0.02 x 10 1 

105.8 

48 Ti 

150 - 1000 

0.08 x 10 10 

265 - 108 

56 Fe 

50 - 1000 

0.2 x 10 10 

832 - 150 

84 Kr 

383 


403 

131 Xe 

228 


1204 

181 Ta 

292 - 313 


1827- 1896 

197 Au 

76 - 165 

lxlO 7 

4828 - 3066 

Sequential Field (Fe/H) 

1000 

Various 

150/0.2 

SPE 

30-180 

Various 

1.26-0.21 


http://www.bnl.qov/nnedical/NASA/CAD/Beann Ion Species and Enerqy.asp 
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Range (mm) 


LET-Range Distribution for Water 



LET (keV/fim) 


http://www.bnl.qov/medical/NASA/CAD/LET-Ranqe.pdf 
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HW1: Mono-energetic Beam 


The variety of ion beams with various energies are used at the NSRL, which 
are relevant to the dominant ion species in space. Exercise physical, 
radiobiological, and nuclear interaction properties of various ion beams for 
materials. 

> Calculate values of LET and range for a water, and compare them to the Range vs. LET 
graph (slide 28). Generate the same graph for another material. 

> Calculate mean free path, probability to suffer a nuclear interaction after 1 or 5 g/cm 2 of 
a material for specified beams and energies. 

> By applying cellular track model of Katz, calculate cell survival probability, translocation 
frequency, HPRT mutation frequency, or HG tumor prevalence frequency from exposure 
to a selected beam. 

> Evaluate radial dose of ionization and excitation for the tissue equivalent material. 

> Calculate the Poisson distribution of ion hits for a given cell size and dose. 

> Evaluate the frequency distributions of energy imparted per DNA target from the direct 
interactions of primary ions (ion events) with the target and the 100 keV electrons (5-ray 
events) produced about an ion's path. 
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HW2: Transport in Thick Target 

The beam delivered to radiobiology samples at NSRL can have many 
components in addition to the nominal beam particles due to the NSRL beam 
transport after beam line shift. Biological effects are determined by both the 
physical beam transport though the targets and the biological effectiveness of 
the mixed charged-particle radiation field. Therefore, biological end points can 
not be described by LET alone. Using the cellular track model of Katz for the 
mixed-radiation fields: 

1) Compare the biological effects as a function of the exposed dose by a given ion 
beam or y-rays after transported through the various depths of tissue equivalent 
material. 

2) Compare the biological Bragg curve of cell death of a given ion beam to the 
physical Bragg curve for a tissue equivalent material. 

> Before the particles reach the Bragg peak region, does the biological Bragg curves follow 

the physical Bragg curve? 

> At the physical Bragg peak location, may not the same peak be observed from the 
biological response curves? 
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Cell Death for Entrance Dose of 1 Gy by 28 Si 300 MeV/u 



> Radiation damage to the DNA is usually the initial event for many radiation induced biological effects observed in 
cells. Most biological end points are usually observed in lightly damaged cells, and the cells heavily damaged and 
unable to replicate will be excluded from analysis for end points. Before reaching the physical Bragg peak, the 
cells are more likely to be lightly damaged by long-range 5-rays. 

> At the Bragg peak region: The particles lose energy sharply and produce 5-rays that have shorter ranges. The cells 
in the physical Bragg peak region are then either heavily damaged when they are directly hit by the charged 
particle or experience less damage by the 5-rays when they are not traversed directly. 

> Cell death curve shows that the severely damaged cells at the Bragg peak are more likely to go through 
reproductive death, the so called "overkill". Not the same peak observed 
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NSRL GERMcode GUI Outputs 


Pv 


NASA SPACE RADIA 



PN [gROCRAM 


The GCR Event-based Risk Model (GERM) code provides scientists, who participate in experiments at the NASA Space 
Radiation Laboratory (NSRL) or in data interpretation of such experiments, the ability to model the NSRL beam line, 
shielding of samples and sample holders, and to estimate basic physical and biological outputs of such experiments. It 
is developed by the NASA Space Radiation Program at Johnson Space Center. 

NASA Space Radiation Laboratory Pictures. 





Brookhavan National Laboratory 


NSRL 


Optical bench in NSRL target room 


Developed By NASA 2011 
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Option for Mono Energetic Beam 
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GERMccce vl.l 2011 
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LET 
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LET and Range curves of 28 Si beam on aluminum shielding as a 
function of kinetic energy of the beam. 
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Option for Mono Energetic Beam 


GERMccce vl.l 2011 
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Mono Energetic Output 

Beam Transport Output 


Options 



28 si(600 MeV/u) on Aluminium 



Energy, MeV/u 


m.f.p 

P(x-l) 

P(x=5) 


Nuclear Extinction: The mean free path and the probability of 
nuclear interaction at 1 g/cm : and 5 g/cm 2 of 28 Si beam on 
aluminum shielding as a function of kinetic energy of the beam. 
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Option for Mono Energetic Beam 


GERMccce vl.l 2011 


S3 

<S> 


■ k~l 

& 

LET Range Curve 
Nuclear Extinction 

& Biological Damage 
Probability of Hits 

Depth Dose g* Multiplicity ai 



& 


& 



o 

|«o| 

Mono Beam 



View 

Print 

Save 

Energetic Transport 

03 

Radial Dose 

& Energy Deposition of DNA 

3 ? Bragg Curve 


Data 


Graph 

Models 


Mono Energetic Output 

Beam Transport Output 



Options 
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• C 


28 si(600 MeV/u) on Tissue 
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Radial Dose: The evaluation of radial dose of ionization, excitation and 
total as a function of radial distance from the exposure to 600 MeV/u 28 Si 
beam on the tissue equivalent material. 
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Option for Mono Energetic Beam 



Biological Damage: Cell survival responses as a function of the dose 
using 600 MeV/u 28 Si beam or y-rays on the tissue equivalent material 
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28 si(60Q MeV/u) on Tissue 

Probability of hits (LET approximation ignoring delta-rays) 



DOSE(Gy), Fluence 
(l/micron2), LET 
(keV/micron), Z*2/beta2 

1.00 0.124E+00 

50.39 312.08 

Mean Hits per Cell by Track- 
Core or with Delta-rays (>1 
mGy) 


Probability of Hits: The Poisson distribution of probability of hits from the exposure to 1 Gy 
using 600 MeV/u 28 Si beam in the cell area of 100 pm 2 . Displayed figure is the probability of 
hits per cell in the track core only, ignoring 5-rays. The mean hits per cell for LET 
approximations with and without 5-rays are stated in the text of the output window. 
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Mono Energetic 

LET Range Curve 

Nuclear Extinction 

Radial Dose 

Biological Damage 
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Energy Deposition of DNA i 




Energy Deposition in DNA Volume: The frequency distributions of energy imparted per DNA volume of 
nucleosome from the exposure to 1 Gy using 600 MeV/u 28 Si beam. The nucleosome (160 Base-pairs) 
volume is assumed as a cylinder in the dimension of 10 nm (diameter) x 5 nm (length). 

• the direct interactions of primary ions (ion events) with the target 

• the 100 keV electrons (5-ray events) produced about an ion’s path 
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Beam Transport 

Depth Dose 


< ► C] 



Depth Dose for Fixed Depth: Normalized depth-dose evaluation using 600 MeV/u 28 Si 
beam on water for primaries, fragments, and total at 10 depths of water in the NSRL beam 
line. The last column shows the fluence-based average LET at each depths. 4 
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DPri = 1.389 

DHI Frag = 0.694 

DTot = 2.083 

<LET> = 103.883 keV/pm 


28 si(600 MeV/u) on Water 
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Bragg Curve: Normalized dose of 600 MeV/u 28 Si beam for primaries, fragments, and 
the total as a function of depth of water in the NSRL beam line. Variations of the 
normalized doses and the fluence-based average LET with water depths are displayed in 
the inset, as cursor moves along the graph. 
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Charge Distribution for depths 


28 si(600 MeV/u) on Water 



■ 4.325 g/cm 2 

• 8.649 g/cm 2 

* 12.972 g/cm 2 



Charge Distribution: Cumulative spectrum of fragments from exposure to 600 MeV/u 28 Si 
beam transported through the water. 

Various number of sets are available in the drop-down menu for the shielding depths by 
selecting “For Fixed Depth”, “For Bragg Curve”, or “For Biological Model” . 42 
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28 si(600 MeV/u) on Water 



Multiplicity of Events: Multiplicity of a-particles in heavy ion fragmentation event at the depth of ~21 g/cm 
of water from exposure to 600 MeV/u 28 Si beam. Drop-down menus of “Depth” and “Ion” are available for 
the selection of depth and the particle type of neutron, proton, or a. Also, energy distribution of heavy ion 
event is displayed by selecting “HI Event” from the drop-down menu of “Ion”. 
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28 Si(600 MeV/u) on Water 
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Multiplicity of Events: Downgraded energy distribution of particles produced from heavy 
ion events at the depth of ~21 g/cm 2 of water from exposure to 600 MeV/u 28 Si beam is 
displayed by selecting “HI Event” from the drop-down menu of “Ion”. 
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Depth I 



28 si(600 MeV/u) on Tissue 



y-Rays 
HZE Nuclei 


Biological Damage: Harderian gland tumor prevalence curve at the depth of ~6 
g/cm 2 of tissue equivalent material as a function of the dose using 600 MeV/u 28 Si 
beam or y-rays. 
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NSRL Facility 

RacfcabonEP91-94.pdf 

Stopping Power and Range Database 

NASATP_LET.pdf 

QMSFRG Model 

RacfcabonEP101.pdf 

Isotopic dependence of GCR 

RacfcabonM.pdf 

Alpha-Nucleus Cross Section 

tp3285_HEdatabase.pdf 

Alpha-Cluster Cross Section 

PRC_cudnotta.pdf 

Radial Dose Model 

tp3497_199S.pdf 

Heavy Particles in Therapy 

Katz_Heavy_Therapy.pdf 

Heavy Ion Mutation I 

NasaTP_Mut.pdf 

Heavy Ion Mutation n 

URB_Mutabon.pdf 

Cellular Track Model 

Katz_NASATP_ 199 1 .pdf 

Frequency Distributions in DNA from HZE Particles 


Event-Based Risk Models 

Event-Based Risk Modeis.pdf 

Haderian Gland tumor Dose Response Models 

MUT10901.pdf 


Opening of “File” tab in the main toolbar. Detailed information about GERMcode 
can be accessed from the menu listed in the left panel of the window. From the 
menu of “Papers”, the published papers can be accessed for the reference to 
GERMcode. 


46 



